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Abstract 

Precise determinations of the strong coupling constant in liadronic collisions demand 
sets of parton densities spanning a sufficiently wide range of values for the QCD scale 
parameter A. For such applications, we supplement the recent GRV parametrization 
by five sets of dynamically generated parton distributions for A-values corresponding to 
as{M'^) between 0.104 and 0.122. Present DIS data, including recent small-x HERA 
results, do not appreciably discriminate between these sets. 
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Parton distributions are essential for connecting observables measured in hard hadronic 
collisions to the underlying partonic interactions. Parametrizations of these quark and 
gluon distributions always include a specific value of the QCD scale parameter A via the 
appearance of this quantity in the renormalization group evolution equations. This sim- 
ple fact leads to a question of consistency when such parametrizations are employed 
in extractions of the strong coupling constant as{Q^) from data on, e.g., jet produc- 
tion rates at the DESY ep collider HERA p| or on l^+jet production at the Fermilab 
TEVATRON Hence precise analyses of this type demand sets of parton densities 
spanning a sufficiently wide range of values for A. Recent parametrizations of CTEQ 
MRS [^, ^ and GRV [0, ^, however, have been provided only for some 'central' values, 
A^ = 200 .. . 250 MeV, of the NLO (MS) scale with four active quarks flavours. The 
only exception is CTEQ2 [Q, where one additional large-A set has been included. 

In this paper, we present five parton density sets which cover, together with the 
recent GRV parametrization [§], the range of A^ from 150 to 400 MeV in steps of 50 
MeV, corresponding to as{M'^) between 0.104 and 0.122. For the derivation of these 
sets we follow the approach developed in 0, in which the quark and gluon densities 
are generated from valence-like inputs at a very low resolution scale = fi"^. Within 
their uncertainties due to the precise choice of fi^ 0, those 'dynamical' distributions 
quantitatively anticipated the strong small-x rise of the proton structure function 
recently observed at HERA |jlO|, |ll|]. Concerning the present task, it should be noted that 



this framework provides a close connection between the experimentally tightly constrained 
quark densities and the much less well-known gluon distribution. This connection leads 
to parametrizations which are as similar as possible for different values of A and thus 
are very well suited for studying the quantitative importance of the consistency question 
raised above. 



Technically, the present paper closely follows the recent GRV update [^] wherever 
possible. Our analysis is carried out in the next to leading order (NLO) perturbative QCD 
framework [Q. We employ the MS renormalization and mass factorization schemes with 
a 'running' number / of active flavours in the QCD beta function and a flxed number 
of three light partonic quark species present in the Q^-evolution ||l3l. The NLO relation 
connecting as(Q^) and Aj^ is chosen as in eq. (10) of p], and the heavy quark thresholds 
in the beta function are taken to be Qh = Trih = 1.5 (4.5) GeV for the charm (bottom) 
case, respectively. The values of Aj^, connected by requiring the continuity of as at 
these thresholds, are listed in Table 1 below. The charm contribution to deep inelastic 
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(DIS) structure functions is approximated by the leading order (LO) massive photon-gluon 
fusion ('Bethe-Heitler') process jl^, since a computationally sufficiently fast version of 
the next order results of []T3[, although already announced |T^, is not yet available. Here, 



we use the LO GRV gluon density with A^q = 200 MeV, and the renormalization and 



factorization scale is fixed at 4m^, irrespective of the value of |13]. For calculating 
fixed-target and ISR inclusive prompt photon production cross sections the complete NLO 
program of [l^ is utilized, with the photonic fragmentation functions taken from . The 



small heavy quark contribution to this process is neglected. 

The construction of our distributions proceeds in two steps: first we fix the non-singlet 
quark combinations = u — u, = d — d, and A = d — u hj adopting and suitably 
modifying the corresponding results of a recent global fit @, then the free parameters 
of the valence-like sea quark and gluon input densities are determined from data on the 
proton and deuteron structure functions ''^ ||19|, [2^, |2l|] and on direct photon production 



in pp collisions |2^, At a reference scale of Qq = A GeV^, the non-singlet densities 



f = Uy, dy, and A are parametrized as 

xf{x,Ql) = AfX^'fil-x/fil + tfy^ + ^ffx) . (1) 

For A^ = 200 MeV they have been directly adopted from the MRS (A) global fit in 
P). Taking into account that the factorization scheme of slightly differs from the one 
employed in 0, no readjustment for the small A-difference is necessary. For the new sets, 
corresponding distributions are derived in the following way: the contributions of u^, dy, 
and A to F2 are calculated for A^^ = 200 MeV in the kinematical region covered by 
the fixed-target data employed in the MRS (A) analysis, and the new inputs for A^^ 7^ 
200 MeV are fitted to these results. For the present purpose, this procedure constitutes 
a sufficiently accurate substitute for a direct fit to data. The resulting parameters for the 
input (1) are presented in Table 1. 

We now turn to the gluon and sea quark input distributions. Here we first have to fix 
the low resolution scales /i^(A) at which the valence-like boundary conditions are imposed. 
In a LO analysis, /x^, Q^, and A enter the parton density evolution and the DIS structure 
function F2 only in the combination as{^Ji'^) / otsiQ"^)- Therefore, if the quark and gluon 
densities are forced to be quite similar at a certain 'central' scale Q\ by data, then also 
the inputs at //^ have to be nearly the same irrespective of A, provided /i^(A) is chosen by 
keeping as{^^^) / (^s{Q\) constant when A is changed. In the NLO the situation is slightly 
complicated by the additional appearance of A via asiQ"^) in the evolution equations and 
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Table 1: The NLO (MS) QCD scales for / active quarks flavours, A^g^, and the indepen- 
dent non-singlet input parameters for eq. (1). The valence quark normalizations, and 
are flxed by the quark number sum rules. 



the expressions for ¥2- Such effects, however, are already incorporated in the A-dependent 
valence distributions derived above, suggesting that some valence moment should be kept 
flxed in this case. In fact, adopting the NLO valence quark momentum fraction < x >„ 
at \}} from <x>.u {fi'^) = J^dxx{uy + ^ 0.58, leads to very similar small-x 

predictions in the HERA kinematic regime for the different choices of A as illustrated in 
Figure 6 below. The resulting scales /i^(A) are listed in Table 2. Their precise values 
are, of course, subject to a similar uncertainty as the value for A^^ = 200 MeV, 
/i^ = 0.34 GeV^ I]. 

At these low resolution scales /i^, the gluon and antiquark distributions are paramet- 
rized as 

xg{x,fi'^) = Agx'^^il-xf^ 
x{u + d){x,fx^) = A^x°i{l-xf^ (2) 
xs{x,fi'^) = xs{x,fi'^) = . 



3 



A^/MeV 


150 


250 


300 


350 


400 


/i / Vjrt! V 


n 99^; 
u. zzo 




n fi9 


n 77 


u.yo 


<X>v 


0.579 


0.582 


0.579 


0.580 


0.582 


ag 


1.9 


1.9 


1.9 


1.9 


1.9 




4.0 


4.0 


4.0 


4.0 


4.0 


2 <x>^_,_j 


0.140 


0.145 


0.147 


0.150 


0.153 




0.35 


0.33 


0.30 


0.29 


0.27 




7.25 


6.84 


6.62 


6.42 


6.20 



Table 2: The input scales /i^ and the independent sea and gluon parameters for the 
valence-hke input in eq. (2). The antiquark normahzation is expressed via the total sea 
momentum fraction 2 <x >u+d) the gluon normalization is fixed by the momentum sum 
rules. 



Our ansatz for the light sea input is slightly simpler than the one in since we do not 
impose constraints from Drell-Yan lepton pair production in pp interactions on the very 
small sea quark densities at large x. The resulting difference of the present sets to the 
NLO GRV parametrization P] is immaterial, except for very special observables. Due to 
the energy-momentum sum rule, the ansatz (2) consists of five independent parameters, 
which are inferred from DIS structure function and direct photon production data. The 
sea parameters are mainly determined by the Ff''^ results of BCDMS |jl9| (normalized 
down by 2.5%), SLAG and NMC 0. This data is fitted in the region where these 
structure functions are sensitive to the sea and gluon distributions, x < 0.3, and where 
higher- twist contributions are expected to be small, > 5 GeV^. The large-x gluon 
density is mainly constrained by the inclusive prompt photon production data of WA70 
|2^ and R806 [^. Here, we fix the factorization scale at = 0.5 where pr denotes 



the photon transverse momentum, but we allow the renormalization scale /i^j to vary with 
A. It turns out that adopting the gluon shape parameters from the NLO (A.^ = 200 

iVi o 

MeV) analysis of leads to a very good description of the data also for the other A-values 
considered in this paper. We have checked that this agreement can be only marginally 
improved by fitting Ug and (3g, thus their values are kept fixed. The sea and gluon input 
parameters for eq. (2) are listed in Table 2. 

A detailed comparison of our results to the data used for fixing the valence-like inputs 
(2) is presented in figs. 1 and 2. We have selected the three sets of distributions derived 
for A = A^^ / MeV = 150, 250, and 350 for the figures. In fig. 1 we compare our fitted 
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results with the BCDMS, SLAG, and NMC fixed-target data |19, 20, 21| for at a; < 0.3 



and >h GeV^, the agreement with the corresponding F2 results is very similar. In fig. 
2, the pp prompt photon data of WA70 (fixed-target) and R806 (ISR) is compared with 
our distributions via the 'default quantity', {crexp — o'th)/crth, with (Xth given by the A = 250 
set. All these comparisons exhibit an excellent agreement with the data, irrespective of 
the value of A. 

The valence-like gluon and sea input and the valence quark densities at = fi'^ are 
shown in fig. 3. The similarity especially between the quark inputs for different A is 
striking; it should be kept in mind, however, that /i^ increases by a factor of three when 
Aa is raised from 150 to 350 MeV. The evolution of these inputs to scales of interest 
for as determinations is illustrated in fig. 4 and fig. 5 for the gluon distribution and 
the quark densities {u and m), respectively. Finally, fig. 6 shows the dynamical small-x 
predictions for F2 in comparison with recent data from the ZEUS and HI collaborations 



|[T0| , 11]. Also shown are the final small-x small-Q^ results of E665 [^. Within the present 
experimental accuracy all sets derived in this paper are in quantitative agreement with 
the HERA results; the for the 149 HERA F^ data points lies between 100 and 110, 
when the systematic and statistical errors are added quadratically and the normalization 
uncertainties are neglected. Future data may impose substantial constraints, especially 
at small Q^,Q^ ^3 GeV^ and very small x where the A-effects are most pronounced. 

To conclude, we have derived five A-variants of the NLO 'dynamical' GRV parton 
densities, together with the standard set spanning the range from 0.104 to 0.122 for 
as{Mz), for use in 05 determinations from hadronic collisions. The step size in as{Mz) 
is about 0.03, which should be sufficiently small for all applications. It remains to be seen 
whether experimental analyses will soon reach an accuracy that renders the inclusion of 
A-differences in the parton distributions mandatory. The Fortran code for calculating 
these parton densities is available via electronic mail from avogt@x4u2.desy.de. 
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Figure Captions 

Fig. 1 The proton structure function Ff' for our parton density sets with A = / MeV 



= 150, 250, and 350 compared to the data at medium x and large ]TP| , |2C 
employed in the fit of our sea and gluon inputs (2). 



Fig. 2 Comparison of the direct photon production data |p2| , P3[[ used for constraining 
the gluon inputs in (2) with the results of three representative parton density sets. 
The factorization scale of the calculations is fip = Pt/2. Data and curves have been 
normalized to the A = 250 result, xt = 2pt/\^, where ^/s is the CMS energy. 

Fig. 3 The valence-hke input distributions xf (/ = m^,, d^j, u, d, g) at the small scales iJ? 
given in Table 2. The strange sea s = s vanishes at = ^ . 



The results at = 5 and 10^ GeV^ have been multiplied by the numbers indicated. 



Fig. 4 The Q^-evolution of our gluon distributions for A^^ = 150, 250, and 350 MeV. 



Fig. 5 The same as Fig. 4, but for the u and u quark densities. For clarity, the antiquark 
distributions are only shown where they substantially differ form the quarks. 

Fig. 6 Comparison of our dynamical leading-twist small-x predictions for with recent 
ZEUS [|ig and HI ^ HERA data and low-Q^ E665 [|3 data. Note that these 
results have not been used in our fits. 
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Fig. 4 
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Fig. 5 
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Fig. 6 
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